Critical dynamics of a spin-5/2 2D isotropic antiferromagnet 
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We report a neutron scattering study of the dynamic spin correlations in Rb2MnF4, a two- 
dimensional spin-5/2 antiferromagnet. By tuning an external magnetic field to the value for the 
spin-flop line, we reduce the effective spin anisotropy to essentially zero, thereby obtaining a nearly 
ideal two-dimensional isotropic antiferromagnet. From the shape of the quasielastic peak as a 
function of temperature, we demonstrate dynamic scaling for this system and find a value for the 
dynamical exponent z. We compare these results to theoretical predictions for the dynamic behavior 
of the two-dimensional Heisenberg model, in which deviations from z = 1 provide a measure of the 
corrections to scaling. 
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Recently, interest in the two-dimensional (2D) square- 
lattice Heisenberg antiferromagnet has intensified due in 
large part to the discovery of high-temperature super- 
conductivity in the doped lamellar cuprates and the sub- 
sequent realization of the near-ideal 2D Heisenberg na- 
ture of their parent compounds £1 The nearest- neighbor 
Heisenberg model is defined as: 

H = J S i • S i C 1 ) 

<i,j> 

where J is the nearest neighbor coupling which is pos- 
itive for an antiferromagnet. Classically, Sj is a three 
component vector of magnitude \J S(S + 1) representing 
the spin at site i, while quantum mechanically Si is the 
quantum spin operator. 

As a result of a symbiotic interplay among theory, 
simulation, and experiment, great progress in under- 
standing the instantaneous spin correlations of the 2D 
Heisenberg antiferromagnet has been made in recent 
years. Chakravarty, Halperin and Nelson (CHN)la de- 
veloped an effective field theory from which an exact 
low-temperature expression for the instantaneous corre- 
lation length, £, has been foundH. While this expression 
agrees closely with experiments on spin- 1/2 Heisenberg 
systemsQ, measurements on systems with S > 112 display 
strong deviations from the predicted behavioioO. Subse- 
quent workErQ has pointed towards a broad crossover from 
classical behavior at high temperature to the renormal- 
ized classical regime where field theory is valid. Recently, 
HasenfratzEl incorporated cutoff effects into the field the- 
ory formalism to describe the behavior in this crossover 
region. 

The dynamics of the 2D Heisenberg antiferromag- 
net has likewise been the subject of detailed theoreti- 
cal work. Tyc, Halperin and Chakravarty (THC)u com- 
bined rcnoxmalization group analysis and dynamic scal- 
ing theoryli3 with simulations of the classical lattice ro- 
tor model to predict a form for the dynamic struc- 



ture factor. Classical molecular dynamicsEj and quan- 
tum Monte Carlo simulations^ have lent credence to 
their predictions; however, due to a lack of suitable sys- 
tems, comparatively few experimental studies of dynam- 
ics in 2D Heisenberg antiferromagnets have been per- 
formed. Some of the most ideal 2D Heisenberg sys- 
tems (La2CuC>4 and Si^CuC^C^) have a very large in- 
tersite coupling J (J w 1500-ftT), making quantitative 
results using conventional neutron scattering techniques 
difficult to obtain. Consequently, previous experiments 
have not resolved the quasielastic scattering from the 
long-wavelength spin- wave excitationsO'tJ. In this com- 
munication, we present a neutron scattering study of 
Rb2MnF4, a quasi- two-dimensional spin-5/2 system with 
an effective spin anisotropy that can be tuned to zero us- 
ing an external magnetic field. Our results provide a de- 
tailed characterization of the dynamic structure factor in 
the quasielastic region. Previous studies&EEl indicate that 
this system behaves like a nearly ideal 2D Heisenberg 
antiferromagnet. Accordingly, we compare our findings 
with the current theoretical understanding of 2D Heisen- 
berg critical dynamics. n 

Following a strategy introduced in our previous workO, 
we exploit the presence of a bicritical point in the field- 
temperature phase diagram of Rb2MnF4 to make possi- 
ble a study of the dynamic spin correlations of a near- 
ideal Heisenberg system over a large range of correlation 
lengths. Rb2MnF4 has the tetragonal K2NiF^ crystal 
structure with an in-plane lattice constant of a = 4.215 A 
and an out-of-plane lattice constant of 13.77 A. The large 
ratio of the out-of-planc to the in-plane lattice constant 
combines with the frustration due to the body-centered 
stacking to make it a nearly two-dimensional magnetic 
system, with an interplane coupling of less than 10 -4 J. 

At zero field, Rb2MnF4 is, a weakly Ising antiferro- 
magnet with J = 0.63 meVO. This interaction energy 
J is more than two orders of magnitude smaller than 
that of the lamellar copper oxide Heisenberg systems, 
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thus making the energy scale of the dynamics much more 
accessible for neutron scattering studies. The principal 
spin anisotropy is a.jnagnetic dipole interaction, with 
g/isHA = 0.032 mcVLJ along the c axis (perpendicular to 
the magnetic plane) . Correspondingly, when a field of ap- 
proximately 5.5 T (depending on temperature) is applied 
parallel to the c axis, the spins flop into the plane. Above 
this spin-flop transition, the system has XY symmetry. 
Precisely along the spin-flop line, and on the extension 
of the line into the paramagnetic phase, the anisotropy is 
effectively zero, so that the system should be in the 2D 
Heisenberg universality class. (See Fig. 1.) 

Experiments were conducted at the NIST Center for 
Neutron Research using NIST's 7 Tesla superconducting 
magnet. We aligned the c axis within 0.5° of the mag- 
netic field to minimize any induced in-plane anisotropies. 
We took field scans at several temperatures to confirm 
the phase diagram and found the line of zero anisotropy 
in Tesla (shown in Fig. 1) to be approximately: H = 
V28.09 + 0.23T where T is the temperature in Kelvin. 
This is in accordance with the form given by Cowley et 
alB. 

Studies of the quasielastic scattering were performed 
with the thermal neutron triple-axis spectrometer BT9 
and the cold neutron spectrometer SPINS. At BT9, we 
used a fixed initial energy of either 13.7 or 14.8 meV 
with a pyrolytic graphite filter before the sample to re- 
move higher harmonics in the incident beam. Collima- 
tions of 40'-27'-Sample-24'-60' were typical, giving an en- 
ergy resolution of 0.8 meV full width at half maximum 
(FWHM). For lower temperatures where higher resolu- 
tion was needed, we used SPINS with a fixed final energy 
of 4 meV and collimations of guide-20'-S-20'-open, which 
gave a resolution of 0.12 meV FWHM. 

Figure 2 shows the scattered intensity as a function of 
energy at the antiferromagnetic zone center for several 
temperatures. The two-dimensional Heisenberg antifer- 
romagnet, in accordance with the Hohenberg-Mermin- 
Wagner theorem, has no transition to long-range order 
above zero temperature. At non-zero temperature it has 
correlated regions whose characteristic length scale di- 
verges exponentially with inverse temperature. These 
correlated regions have a finite lifetime which translates 
into a non-zero energy width of the quasielastic peak pro- 
duced in the dynamic structure factor. As the temper- 
ature is lowered towards zero, the correlated regions be- 
come progressively more stable, and the energy width of 
the peak decreases. The measurements in Fig. 2 display 
this critical slowing down. 

According to dynamic scaling theory, the functional 
form of the structure factor is independent of tempera- 
ture. The temperature dependence enters only through 
the reduced reciprocal space position k (2D reciprocal 
space distance from the magnetic zone center) and fre- 
quency u>, which are scaled by £, the correlation length, 
and ujq, the characteristic frequency, respectivelytHl: q = 
kt; and v = uj/loq, so that q and v are both dimensionless. 
In addition, the characteristic frequency is predicted to 



scale with the correlation length to a power — z, with 
z = d/2 for a Heisenberg antiferromagnet, where d is the 
spatial dimension. 

In accordance with these predictions, we fit the energy 
scans through the quasielastic peak at (0 1 0) to the 
dynamic structure factor: 



S(k 1 u)=^ 1 S(qMq,iy). 
We took Lorentzian forms for S(q) and &(q, v): 

So 



S(q) = 



$(q,v) 



1 + q 2 
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(/z is an arbitrary constant.) Due 



to the finite resolution, g-dependent contributions were 
needed to reproduce accurately the observed lineshape. 
/j, = 1.7 ± 0.2 gave the best fit at all temperatures, in 
agreement with the values (1.4 and 2.0) found by THC in 
their analyses of the classical lattice rotor model. For the 
fits, £ n was fixed at the values determined in our previous 
studyO, and the temperature dependence of So was found 
to agree well with the results of those measurements. The 
correlation length for the temperature range accessible 
for studying the dynamics varied from 1 < £/a < 60. 

As Fig. 2 demonstrates, fits to this form convolved 
with the experimental resolution are quite good at all 
temperatures. Figure 3 shows the results for the energy 
widths extracted from these fits. Note that data taken 
at SPINS and BT9 agree closely in the overlapping re- 
gion. This agreement gives us confidence that we have 
correctly accounted for the very different experimental 
resolutions in the two measurements. E3 When the tem- 
perature is scaled by JS(S + 1), the temperature de- 
pendence of u>o(T) agrees well with the results of classi- 
cal molecular, dynamics simulations carried out by Wysin 
and Bishopllil. They predict a temperature scaling factor 
of JS 2 , but when scaled by this factor, our data show a 
much stronger temperature dependence than that exhib- 
ited by the simulation data. Normalizing temperature 
by the classical spin stiffness JS(S + 1) has been showncl 
to collapse the instantaneous correlation length data for 
2D quantum Heisenberg antiferromagnets with S > 1 
onto the classical results at high temperatures, and here 
again succeeds in reconciling the spin-5/2 data with cor- 
responding results for the classical system. Thus, as with 
£, the dynamic behavior follows classical scaling at high 
temperature (1 < £/a < 10). 

Similar measurements of the quasielastic energy width 
have been performed by Fulton et a/E3 on KFeF 4 , an- 
other 2D spin-5/2 antiferromagnet. These results agree 
with our data at the highest temperatures, but deviate 
strongly at lower temperatures. We believe that this dis- 
crepancy results from a crossover to Ising criticaLbehav- 
ior in KFeF 4 . Studies of Rb2MnF 4 in zero fielcH show 
that the Ising crossover occurs near 1.2Tjv- KFeF 4 has 
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nearly the same reduced Ising anisotropy as Rj^MnF/t, 
and hence would also be expected to enter a region of 
Ising critical behavior below 1.27V. All but the highest 
temperatures from the KFeF4 study therefore lie below 
the Ising crossover region. 

As mentioned above, dynamic scaling theory pre- 
dicts ojq oc £~i_with z = 1 for the 2D Heisenberg 
antiferromagnetll3. In addition, CHN predict corrections 
to scaling which go as T 1 / 2 : lj = c ^ 1 {j§]r) 1 ^ 2 where c 
is the spin-wave velocity, and p s is the spin stiffness. Fig- 
ure 4a shows a plot of ojq versus l/£; the best fit to the 
simple form ljq oc gives z — 1.35±0.02. This value for 
z is intermediate between the values for the 2D (z — 1) 
and 3D (z = 1.5) Heisenberg antiferromagnets. However, 
as detailed below we believe only a 2D model is relevant 
here. Figure 4b, which shows the product ujq^ versus 
temperature, demonstrates the corrections to scaling if 
z = 1 is assumed. Clearly, corrections are stronger than 
the T 1 / 2 predicted by CHN. Simulations of the classical 
model, as mentioned above, agree with our data for luq, 
yet they claim to see a different temperature dependence 
for the product wo£- This is most likely due their use of 
a form for £(T) that has since been shown to be inaccu- 
rate in this temperature range. Monte Carlo studies on 
a spin- 1/2 systemEj have also indicated that z = 1, but 
with a temperature independent prefactor ip agreement 
with predictions by Arovas and Auerbach.E3 

The data in Figs. 3 and 4, taken at face value, may 
suggest that z > 1 for Rb2MnF4 along the bicritical line, 
or that there is a crossover to some other critical behavior 
with a non-zero phase transition temperature. Explana- 
tions involving a crossover to three-dimoBsianal behavior 
seem unlikely in light of previous studiescJ'Eil at zero field 
showing that Rb2MnF4 behaves as a nearly ideal two di- 
mensional system to very large correlation lengths. Like- 
wise, 2D Ising or 2D XY behavior are precluded by the 
high temperature at which we observe deviations from 
z = 1 behavior, as compared to the scales at which these 
crossovers should occur, as well as by our previous re- 
sults on the staticsEl, which agree very well with theory 
and simulation for the 2D Heisenberg model. 

However, the dynamic scaling near the bicritical point 
could still conceivably differ from that of the ideal 2D 
Heisenberg antiferromagnet. While the universality class 
for static critical behavior is determined solely by the 
symmetry properties and the spatial dimension, the dy- 
namics can also be affected by conservedrquantities and 
the Poisson-bracket relations they satisfytL2l. The bicrit- 
ical region differs from a true isotropic system due to 
the non-zero, conserved uniform magnetization along the 
applied field— direction. Noting this distinction, Dohm 
and JanssenEj performed a renormalization group study 
of bicritical dynamics in 4 — e dimensions. They found 
that dynamic scaling was obeyed, but that the exponent 
for the 3D bicritical point was larger than that for the 
3D Heisenberg model. To explore the possibility that 
we might be seeing a similar effect in our 2D system, 



we measured loq in zero field at temperatures above the 
Ising-Heisenberg crossover. These results, shown in Fig. 
4a, overlap closely with the data taken at the same tem- 
peratures on the bicritical line. This indicates that, for 
this temperature range, the magnetic field itself is not 
measurably affecting the quasielastic width. Clearly, ad- 
ditional theoretical work on 2D bicritical dynamics and 
corrections to dynamic scaling for the 2D Heisenberg an- 
tiferromagnet would greatly elucidate the findings from 
these measurements. 

With these measurements of the dynamic spin corre- 
lations in Rb2MnF4 near the bicritical point, we have 
provided the first experimental study of the quasielas- 
tic behavior in a 2D isotropic antiferromagnet. These 
results are largely consistent with the current theoreti- 
cal understanding of the dynamics of the 2D Heisenberg 
model, but also raise some questions. The shape of the 
dynamic structure factor in the quasielastic region obeys 
a form consistent with dynamic scaling, and the tem- 
perature dependence of the characteristic frequency ujq 
is consistent with the anticipated form £~ 2 , though with 
z larger than the predicted value z = 1. To establish 
whether the difference in z originates in stronger correc- 
tions to scaling than predicted or indicates a distinction 
between ideal 2D Heisenberg dynamic scaling and the 
dynamic behavior near a 2D bicritical point will require 
further theoretical work as well as experimental studies 
of other Heisenberg antiferromagnets. 
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